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ABSTRACT 

We use a Monte Carlo simulation to calculate the spectra of mildly relativistic 
thermal plasmas in pair balance. We use the exact integral expression for the 
electron-positron thermal annihilation spectrum, and provide accurate expressions 
for the Gaunt factors of electron-ion, electron-electron, and electron-positron thermal 
bremsstrahlung in the transrelativistic temperature regime. The particles are assumed 
to be uniformly distributed throughout a sphere, and the pair opacity is self-consistently 
calculated from the energy and angular distribution of scattered photons. The resultant 
photon spectra are compared with the nonrelativistic diffusion treatment of Sunyaev 
and Titarchuk, the bridging formulas of Zdziarski, and the relativistic corrections 
proposed by Titarchuk. We give a corrected formula for the spectral index resulting 
from Comptonization in the low temperature, low optical depth regime, and show 
where the effects of bremsstrahlung production are important in spectral calculations. 

We calculate allowed pair-balanced states of thermal plasmas with no pair escape 
which include bremsstrahlung and internal soft photons. The results are presented in 
the spectral index/compactness plane, and can be directly compared with observations 
of spectra from AGNs and Galactic black hole candidates. By comparing with X-ray 
spectral indices of Seyfert AGNs and compactnesses inferred from X-ray variability 
data, we find that the allowed solutions for pair equilibrium plasma imply that the 
temperatures of Seyfert galaxies are ^ 300 keV. This prediction can be tested with 
more sensitive gamma-ray observations of Seyfert galaxies. We find that if the X-ray 
variability time scale gives an accurate measure of the compactness, pair-dominated 
solutions are inconsistent with the data. 
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1. Introduction 

Comptonization in mildly relativistic plasmas has long been thought to play an important 
role in producing the continuum emission from compact Galactic X-ray sources (Zel'dovich & 
Shakura 1969; lUarionov & Sunyaev 1972; Felten & Rees 1972). It has also been suspected 
that such plasmas reside near the innermost regions of Galactic black holes such as Cygnus X-1 
(Shapiro, Lightman & Eardley 1976; Sunyaev & Titarchuk 1980). In addition to Comptonization, 
pair production and annihilation may also be important in such plasmas. In fact, the emission 
around 1 MeV reported from the direction of the Galactic center (Riegler et al. 1985) and Cygnus 
X-1 (Ling et al. 1987) was interpreted as radiation from pair-dominated plasmas (McKinley 1986; 
Liang & Dermer 1987). Furthermore, observations of transient features with energies near 0.4 MeV 
observed from Galactic black hole candidates (e.g. Gilfanov et al. 1994) could be the signature of 
pair annihilation, although other interpretations are possible (Skibo, Dermer, & Ramaty 1994). 

Evidence for the existence of mildly relativistic thermal plasmas in extragalactic astrophysical 
environments is provided by X-ray and gamma-ray observations of the Seyfert class of Active 
Galactic Nuclei (AGN). Observations at X-ray energies made with HEAO-l and EXOSAT revealed 
a power law spectral shape with a canonical photon spectral index T ~ 1.7 in the 2-10 keV range 
(Rothschild et al. 1983; Turner & Pounds 1989). More recent analyses of Ginga observations of 
Seyfert 1 galaxies, taking into account the effects of reflection by a cold optically-thick medium, 
imply that the mean X-ray spectral index is closer to L ~ 1.9 — 2.0 (Pounds et al. 1990; Nandra 
& Pounds 1994). At gamma-ray energies, observations of Seyfert galaxies made with OSSE on the 
Compton Gamma- Ray Observatory (CGRO) reveal spectra which generally display exponential 
cut ofl"s or spectral softenings with an e-folding energy of ~ 50 keV - several hundred keV (Maisack 
et al. 1993; Johnson et al. 1994; Zdziarski et al. 1995). These observations, in conjunction 
with the rapid temporal X-ray variability observed from these sources (Barr & Mushotzky 1986; 
McHardy 1989; Done & Fabian 1989; Grandi et al. 1992), strongly suggest the presence of mildly 
relativistic thermal plasmas in the cores of Seyfert galaxies. 

For pair-free plasmas at nonrelativistic temperatures and high optical depths, analytic 
formulae for the Comptonized spectra were derived by Sunyaev and Titarchuk (1980, hereafter 
ST80) assuming the presence of a low-energy photon source with luminosity greatly exceeding 
the intrinsic bremsstrahlung luminosity of the plasma. This problem was treated at relativistic 
temperatures using a Monte Carlo method (Pozdnyakov, Sobol & Sunyaev 1977; 1979; 1983; 
Gorecki & Wilczewski 1984). Analytic corrections to the ST80 calculations at mildly relativistic 
temperatures and low optical depths have recently been given by Titarchuk (1994) and verifled 
using the Monte Carlo method of Hua and Titarchuk (1994, hereafter IIT94). The Comptonization 
of bremsstrahlung in mildly relativistic plasmas including pairs was treated extensively by 
Svensson (1982b; 1984) where the conditions for pair equilibrium (i.e., plasmas where the pair 
production rate equals the annihilation rate) were analytically obtained. In a subsequent series 
of papers, Zdziarski (1984; 1985, hereafter Z85; 1986) calculated the emergent spectra from pair 
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balanced plasmas including both bremsstrahlung and synchrotron soft photons using a Monte 
Carlo simulation. Analytic representations of the photon spectra were presented. 

In the present analysis we calculate a grid of solutions for mildly relativistic (0 = kT/nieC^ ^ 2) 
plasmas in pair equilibrium using an iterative Monte Carlo procedure. We present new formulas 
for the thermal bremsstrahlung Gaunt factors at these temperatures and use the exact integral 
expression for thermal annihilation. Moreover, we calculate the pair opacity due to photon- 
photon interactions using the the angle- dependent photon distribution rather than through an 
angle-averaged approximation. We allow for the presence of a soft photon source in addition 
to the intrinsic bremsstrahlung and annihilation radiation of the plasma. This analysis differs 
from past numerical simulations of plasmas in pair balance which considered only bremsstrahlung 
(Zdziarski 1984; McKinley 1986; Liang & Dermer 1987) or thermal cyclo- synchrotron (Zdziarski 
1986; Kusunose & Takahara 1985) soft photon sources. We consider, however, only steady-state 
pair plasmas (see, e.g., Kusunose 1987 for a treatment of the dynamics of thermal pair plasmas). 
In the pair-free case we compare our results with those of ST80, Z85 and IIT94 and specify the 
limits of applicability for the use of these formulae. The pair balanced solutions that we obtain 
are presented in the spectral index/compactness plane where direct comparisons with data from 
high energy observations of Seyferts are made. 

In §2 we discuss the bremsstrahlung and pair annihilation radiation processes intrinsic to 
mildly relativistic thermal plasmas, and in Appendix A we present formulae for the thermal 
bremsstrahlung Gaunt factors valid at mildly relativistic temperatures. In §3 we describe our 
Monte Carlo simulations, and in §4 we give the results for pair-free plasmas and compare them 
with analytic calculations found in the literature. Our results for pair balanced plasmas are given 
in §5 and compared with data from Seyfert galaxies in §6. We summarize our results in §7. 

2. Radiation and Pair Balance in Mildly Relativistic Plasmas 

In a thermal plasma consisting of electrons, positrons and protons at mildly relativistic 
temperatures (0 ^ 2), the dominant processes are Compton scattering, bremsstrahlung, 
photon-photon pair production, and electron-positron annihilation. Higher order processes such 
as photon-particle and particle-particle pair production, double Compton scattering, and three 
photon annihilation are negligible at these temperatures (Svensson 1984) and will not be included 
in this analysis. In addition to the intrinsic radiative processes of the plasma we include, for 
generality, a soft photon source. Such a soft photon source might arise, for example, if the plasma 
is threaded with a magnetic field giving rise to thermal self-absorbed synchrotron emission (e.g., 
Zdziarski 1986). Other sources of soft photons could involve external radiation fields from accretion 
disks. For optically-thick plasmas the spectrum of the emergent radiation will be modified by 
transport effects. Absorption due to the inverse bremsstrahlung process can be important for 
sufficiently low photon energies. In addition, if the plasma has a finite optical depth to Compton 
scattering, then the scattering of these photons with thermal electrons and positrons will play an 
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important role in determining the emergent spectrum. 

The calculation of the photon production due to binary collisions involves integrations of 
the relavent cross sections over the thermal particle momentum distributions. General integral 
expressions for photon reaction rates in relativistic plasmas have been derived by Weaver (1976). 
Formulas for the spectrum of thermal e+e" annihilation radiation have been derived by Svensson 
(1983) and Dermer (1984; for a numerical treatment, see Ramaty & Meszaros 1981). The photon 
production rate per unit volume per unit energy due to thermal pair annihilation is given by 

1^ dlrilr-l)exp(^-^^^aa{jr), (2-1) 

where x = e/nieC^ is the dimensionless photon energy in units of the electron rest mass, K2{^) is 

— * — * 

the modified Bessel function of the second kind, and 7^ = 7+7- (1 — /3+ • /3_ ) is the invariant Lorentz 
factor associated with the relative velocity of the e~^e~ system. The terms n_ and represent 
the volume densities of electrons and positrons, respectively, and (Ta{jr) is the annihilation cross 
section, given by 

(e.g., Jauch & Rohrlich 1976). Here = e^/rrieC^ = 2.82 X 10"^^ cm is the classical radius of the 
electron. 

Relativistic e^p, e^e^ and e~^e~ thermal bremsstrahlung was calculated by Dermer (1986) 
using the lowest order bremsstrahlung cross sections. Here we give a semi-analytic representation 
of the emissivities in the form 

Si2{x; 0) = nin2Q-^^^x-^ exp(-a;/0)5i2(a;; 0), (2-3) 

where we fit various functional forms for the Gaunt factor gi2{x; 0) to the curves obtained by 
numerically integrating the exact expressions given by Dermer (1986). Fits to the Gaunt factors 
in the transrelativistic temperature regime accurate to better than 8% are given in Appendix A, 
along with asymptotic forms for bremsstrahlung soft photon production at nonrelativistic and 
ultrarelativistic temperatures (see references in Svensson 1984). 

For simplicity we consider a pure hydrogen plasma with spherical geometry and uniform 
matter density. We consider only plasmas in pair balance and assume that pair escape is negligible. 
Note that the pair balance condition is nonlinear, insofar as the photon opacity for pair creation 
is supplied by the photons themselves, which are generated in turn by Compton scattering and 
annihilation of the electrons and pairs. The fundamental input parameters of the problem are 
the temperature 0, the radius R, the soft photon luminosity Ls, the proton density rip, and the 
positron density The condition of pair equilibrium is used to iteratively determine and 
the electron density is fixed by the requirement of charge neutrality, n_ = rip + Since the 
reaction rates and luminosities of binary collision processes between any two species with densities 
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Til and 712 scale as nin2R^i it is more convenient to work with the reduced set of parameters 
0, z = n^/np, the Thomson depth to ionization electrons Tp = npaxR, and the soft photon 
compactness is = Ls(Tt / RnieC^, where cry = 6.652 X 10~^^ cm^ is the Thomson cross section. 
The total Thomson optical depth r = (1 + 2z)Tp. In terms of these variables, R and rip enter 
only through the product Tp = npaxR- The luminosity of the emergent radiation scales as R for 
a given value of Tp, whereas the shape of its spectrum is unchanged. This scaling holds whenever 
most of the soft photon energy density is injected at frequencies well above the free-free self 
absorption frequency. This is the case for dimensionless soft photon energies a; ^ 5 X lO~®r0~'^/'*, 
corresponding to the free-free self absorption frequency for an emission region of size R = 3 X 10^ 
cm, which corresponds to 10 Schwarzschild radii for a 10 Mq black hole. 



3. Description of the Monte Carlo Code 



We calculate the cumulative effects of Compton scattering and 77 pair attenuation self- 
consistently using an iterative Monte Carlo approach. When a photon is Compton scattered, 
the final energy and angle through which it scatters are randomly drawn from the distributions 
determined from the Klein- Nishina cross section and the relativistic Maxwellian distribution of the 
electrons. This is done by first performing a Monte Carlo simulation of single Compton scattering 
of a photon with an electron in the thermal plasma, and constructing scattering energy and angle 
arrays. These arrays are used in a second program, described below, which follows the photon 
propagation in the hot plasma. 

In the program that generates the scattering arrays, a given photon with specified initial 
energy x is scattered by an electron chosen with a random direction and energy jnieC^ drawn from 
a relativistic Maxwellian of temperature 0. If the photon scatters through an angle cosine of /Xg 
then the scattered photon will have the energy given by (e.g., Pozdnyakov et al. 1983) 

(3.1) 



— * 

where fii and fif are the cosines of the angles between the initial electron velocity vector /3 and the 
initial and final photon propagation directions, respectively. The angle through which the photon 
scatters, /Xg, is given by 

fis = m + [(1 - - /^/)]'/' cos(,/.^ - (3-2) 

where and (j)f are the azimuthal photon directions. This scattering event is assigned a weight 
given by 

w{x,Xs,fis) = [1 - PfJ-t)-^- (3-3) 
where the polarization-averaged Klein- Nishina cross section for the scattering of a photon with a 
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moving electron is (e.g. Pozdnyakov et al. 1983) 



da 



2 72(1 - /3/x,)2 V X 



2 r 



— + ^ + 2 — - — 
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with 



and 



X* = 7(1 - PfJ'z)x, 



X*s = - /3fJ-f)Xs 



(3-4) 



(3-5) 



(3-6) 



These weights are summed in bins of scattered photon energy, Xs, and scattering angle cosine, /Xg, 
for a grid of incident photon energies. In this way we obtain the probability P[x, fis,Xs)dxsdfis for 
the scattering of a photon with energy x into the angle cosine range dfXs about the angle cos~^ /Xg 
and into the energy range dxs around Xs by an electron in a plasma of temperature 0. We then 
tabulate the cumulative distributions given by 



and 



F^{xs) = dx'^ j ^ dfi'^ P{x, fi'^, a;'J, 

//is 
^ dfJ''sP{x,fJ-'s,Xs). 



(3-7) 



(3-8) 



The energy of the scattered photon, Xs, and the scattering angle, /Xg, can be obtained from 
two uniform random numbers, ^, 77 G [0, 1] by inverting the cumulative probability distribution 
functions: 



Xg 
fJ-s 



(3-9) 
(3-10) 



For a given simulation we first specify the input parameter set 

Pi = {0, Tp, z, h}. (3-11) 

Photons are then injected into the plasma sphere with random direction and specified spatial 
distribution. We assume a uniform distribution except when we make comparisons with the results 
of Sunyaev and Titarchuck (1980), in which case we use the distribution they employed to obtain 
analytic solutions in the nonrelativistic diffusion regime. This distribution is given by 

where r is the distance from the center of the sphere. The energy of the injected photons is 
drawn from the bremsstrahlung, annihilation and soft photon emissivities. For the soft photons 
we assume either a black body spectrum with temperature kT = 0.1 keV or a mono-energetic 
distribution at energy 0.33 keV. The relative abundances of bremsstrahlung, annihilation and soft 
photons are governed by the input parameter set Pj. 
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Although bremsstrahlung absorption is generally unimportant, we determine the relative 
probability that a photon is scattered, absorbed by the inverse bremsstrahlung process, or escapes 
from the system. If the photon is absorbed we exit the loop and choose a new photon. If the 
photon is scattered, a new photon energy and direction are obtained and the process is repeated 
until either the photon escapes or is lost by free-free absorption. Those that escape are binned in 
energy, building up the emergent spectrum. At each step the energy and direction of the photon 
and the radial location of the scattering event are recorded. From this data we obtain the photon 
density n^'^\r, ft, x). 

In the first iteration 77 pair production is neglected. The simulation is then repeated with 
the inclusion of an approximation to the energy- dependent radially-averaged 77 opacity, a^-^[x), 
obtained from the photon density distribution n^'^\r, ft, x) established in the first iteration. First, 
the optical depth t^'^[x) along the radius of the sphere for a photon of energy x is obtained by 
numerically evaluating the integral 



T 

77 



(1) 



(a;) = ^ dr^dnj^^ dx' n^°\r,n,x'){l - fi)a^^{x,x',fi), (3-13) 



where a^^[x, x', /x) is the total cross section for the process 77 e~^e and is given by (e.g. Jauch 
& Rohrlich 1976) 



,{x,x',fi)=^Trrl{l-/3') 



(3-/3^)ln(i±|)-2/3(2-/32) 



(3-14) 



Here /3 is the velocity of either the electron or positron in the center of mass and is related to x, 
x' and /X by the relation 

We evaluate the integral (3-13) by the Monte Carlo technique and then set 

aW(a;) = tI^]\x)/R. (3-16) 

We then repeat the entire process, this time including the effects of 77 attenuation through 
a^^[x). This results in a new spectrum and photon density distribution n^^\r, fl, x') from which 
we obtain a\~l[x) and so forth. For pair-poor plasmas [z <C 1) the first iteration is always sufficient 
for establishing equilibrium. For pair-dominated plasmas [z ^ 1) we find that pair balance is 
usually obtained at the second or third iteration. 



4. Pair-Free Solutions: Comparisons With Analytic Results 



In this section we give results for the case of a pair-free plasma [z = 0) where the soft photon 
luminosity greatly exceeds that of bremsstrahlung. Thus r = Tp. We inject soft photons with 
energy xnieC^ = 0.33 keV, a random propagation direction, and spatial distribution given by 
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equation (3-12) into a sphere of radius iZ = 3 X 10^ cm. The choice of this spatial distribution is 
made in order to compare with the analytic results of ST80. In addition, we compare our results 
with the analytic bridging formulae of Z85 and the relativistic corrections to the ST80 calculations 
given by HT94. We set z = 0, iZ = 3 X 10^ cm and Ig ^ Ijj^ where the bremsstrahlung soft 
photon compactness is given by 

K, 4.3a/r20^/2(l + 1.20 + 1.80^-^'* + 1.30^ + 1.20^ - 1.50^-^) (4-1) 

for < 1 (Svensson 1982) and aj is the fine- structure constant. This reduces the parameter set 
Pi to {0,r}. 

In Figures la-Id we compare the results of our Monte Carlo simulations for ^-function 
soft photon injection with the analytic results of ST80, Z85 and IIT94 for the teperatures 
= {0.1, 0.2, 0.5, 1.0} and Thomson depths r = {0.1, 0.5, 1.0, 5.0}. All curves are normalized 
such that the integral over energy is unity. As can be seen, in the high optical depth, low 
temperature regime of thermal Comptonization, all analytic fits are in good agreement with the 
Monte Carlo simulation. This is expected because the analytic fits are designed to approach the 
Sunyaev-Titarchuk result derived in the large optical depth diffusion limit. As r becomes much 
less than unity, the ST80 results break down whereas the IIT94 analytic fit still provides a good 
representation, because they modify the escape probability used in the Sunyaev-Titarchuk limit 
to apply to the low optical depth case. The fit of Z85 is poor in the low temperature, low r case 
because he is only considering parameters leading to photon spectral indices softer than T = 1.4, 
which requires large optical depths in the low temperature limit. At temperature = 0.2 (Fig. 
lb), both Z85 and IIT94 are in good agreement with the Monte Carlo results at large optical 
depths, although Z85 is again discrepant at low optical depths. At still higher temperatures shown 
in Fig. Ic and Id, serious discrepancies appear between the analytic fits and the simulation results. 

We, therefore, recommend the use of the model of Sunyaev and Titarchuk (1980) for r ^ 2 
and ^ 0.2. The corrections of Titarchuk (1994) (see also Hua & Titarchuk 1994) extend the 
validity of this model to r ^ 0.2 and <^ 0.3. If the photon spectral index is in the range 
1.4 ^ r ^ 2.5, which occurs for small r and large or, alternatively, small and large r, then 
the representation of Zdziarski (1985) is in agreement with our simulations. We advise caution 
in the use of these formulas when r ^ 1 and 0^1 because at these temperatures and optical 
depths, the solutions may not be in pair balance and bremsstrahlung contributes significantly, a 
situation we consider in more detail in §5. 

In Figure 2 we show the variation of the 2-10 keV photon number spectral index, r(2-10 keV), 
as a function of r for the temperatures = {0.1, 0.2, 0.5, 1.0}. We again see good agreement in 
the large optical depth limit, but large discrepancies are found in the low optical depths for the 
analytic fits that are restricted to the diffusion limit (ST80) or to spectra that are relatively hard 
(Z85). The results of IIT94 are in good agreement except at high temperatures. In the low optical 
depth limit of the = 1 curve, however, the Monte Carlo simulations are very sensitive to the 
injection energy of the photons because the amplification factor is large and the Compton-scattered 
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Fig. la. — The emergent photon spectra resulting from Monte Carlo simulations of the 
Comptonization of monochromatic soft photons with energy 0.33 keV (histograms) compared with 
the analytic approximations of Sunyaev & Titarchuk (1980; solid curves); Zdziarski (1985; dotted 
curves), and Hua & Titarchuk (1994; dashed curves). The dimensionless electron temperature 
= 0.1, and spectra for different optical depths of electrons are shown by the labels in the 
separate panels. 




Fig. lb. — Same as for Fig. la, but with = 0.2. 
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Fig. Ic. — Same as for Fig. la, but with 



= 0.5. 




Fig. Id. — Same as for Fig. la, but with = 1.0. 
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Fig. 2. — Monte Carlo calculations of the 2-10 keV photon spectral index resulting from thermal 
Comptonization of monochromatic photons with energy 0.33 keV (heavy solid curves) are compared 
with the analytic approximations of Sunyaev & Titarchuk (1980; solid curves); Zdziarski (1985; 
dotted curves), and Hua & Titarchuk (1994; dashed curves). Also shown by the dot-dashed curves 
are the low temperature, low optical depth approximations to the spectral index given by equations 
(4-2) with equation (4-5) replacing equation (4.3). 
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Fig. 3. — The variation of the 2-10 keV photon spectral index with the temperature of the soft 
photon blackbody for r=0.01, 0.1 and 1.0. 



- 15 - 



spectrum mirrors the effects of the individual scattering events (see Pozdnyakov et al. 1983). 
This effect is shown clearly in Figure 3, where we plot values of r(2-10 keV) as a function of 
the temperature for the blackbody soft photon distribution scattered by thermal electrons with 
= 0.5. When r = 0.01, we see variations of T by nearly 2 units. 

In the low temperature, low optical depth regime, the spectral index is harder than expected 
using arguments by Zel'dovich (see Pozdnyakov, Sobol', & Sunyaev 1977; Zdziarski 1986; Dermer, 
Liang, & Canfield 1991), which imply that the photon index 

r = 1 - In P(r)/ In A. (4-2) 

Here the energy amplification factor 

^(0) = 1 + ^t^/w/,?^ = 1 + 40 + 1602 (4-3) 

iV2(l/0j 

is the ratio of the mean photon energy after scattering to its initial energy in the limit a; ^ and 

Pir) = 1 - ^[{2r' - 1) + e-2-(2r + 1)] (4-4) 

is the mean scattering probability in a uniform sphere (see Osterbrock 1974). 

We find that the formation of the spectrum in the low temperature, low optical depth regime 
is controlled by the spread of the photon number distribution rather than through the average 
increase in energy of a photon after scattering. Approximating a nonrelativistic thermal electron 
distribution by a monoenergetic distribution with speed /3 = (30)^^^, one finds that with each 
scattering, the photon distribution spreads by a factor ~ 1 + (30)^/^. Using analogous arguments 
to those given by Zel'dovich, this implies a spectral index in this regime given by equation (4-2), 
but with the amplification factor A replaced by a new factor A' ~ 1 + (30)^/^, representing the 
spread of the photon distribution function. We can generalize the spectral index approximation 
(4-2) to arbitrary temperatures by replacing A[Q) with the new factor 

A*(0) ^ 1 + (30)^/2 ^40 + 1602. (4-5) 

Good agreement with the numerical results is obtained in the low temperature, low optical depth 
regime with this replacement. Figure 2 shows a comparison of the r ^ approximation using 
equations (4-2) with equation (4-5) replacing equation (4-3). 

5. Pair Balanced Solutions 

In this section we give results for the pair-balanced solutions obtained with our Monte Carlo 
simulations. In addition to the intrinsic bremsstrahlung and annihilation radiation of the plasma, 
we inject soft photons with a blackbody energy distribution at temperature kT = 0.1 keV and 
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soft-photon compactness Ig. We vary the parameter set P/ = {0, Tp, z,ls\ by first fixing 0, Tp 
and Is and then adjusting z until pair balance is achieved. There are either zero, one or two 
pair-balanced states for a given 0, Tp and Ig (Svensson 1984; Zdziarski 1984, 1985). A low z 
solution occurs when the photon density at energies above threshold for pair production is low 
enough such that the small amount of pair production that ensues can be balanced by the steady 
state annihilation of a relatively small amount of pairs. As the pair fraction z is increased, the 
annihilation rate will at first exceed the production rate until the scattering opacity due to the 
pairs becomes comparable to that of Tp. At this point the spectrum will harden as more low 
energy photons get Comptonized to higher energies. The pair production increases in a nonlinear 
fashion until balance is again obtained. This is the high z, pair dominated regime. 

In Figure 4 we show the results for Ig = 0. The compactness of the emergent radiation Ih 
is plotted as a function of 0, where each curve represents a different value of Tp. We compare 
our results with the analytic calculations of Svensson (1984) and the Monte Carlo simulations of 
Zdziarski (1985) for the case Tp = 1. Factor-of-two deviations between our numerical results and 
the analytic treatment can be attributed to the uncertainty in deriving the number of scatterings 
required for a soft photon to enter the Wien regime and pair produce. We find good agreement 
between our results and the calculations of Zdziarski in the regime where z <C 1, which occurs 
on the lower branches of the curves in Figure 4. Small discrepancies between the numerical 
simulations in the upper pair-dominated portion of the curves are due to an improved expression 
for the e~^e~ bremsstrahlung emissivity and a different method for calculating t^^. In this work we 
used the full angular distribution of the photons to calculate the pair opacity, whereas the photon 
distribution was assumed to be isotropic in the analysis of Zdziarski (1984). 

Figure 5 shows the effects of increasing the soft photon compactness Is on the allowed states 
of a thermal plasma in pair balance. Here we plot the hard compactness Ih as a function of 
for different values of Tp and Is- The outermost curves in the panel correspond to the Ig = 
solutions of Figure 4. The effect of increasing the soft photon compactness is to increase the 
hard compactness and to decrease the maximum temperature Qmaxi^s^Tp) permitted for a steady 
pair-balanced plasma. 

In Figure 6 we show examples of spectra resulting from our simulations. These spectra 
correspond to plasmas in pair-balance with no external soft photon injection [Is = 0) and the 
parameter values = {O.2, 1.0} and rp={0.1, 1.0}. For each set of parameters we show the spectra 
for the two states of the plasma permitted by the condition of pair-balance. As can be seen, in 
the \ow-z, low optical depth (rp = 0.1) states, the plasma simply emits a thermal bremsstrahlung 
spectrum. The dotted curves for the z = cases correspond to pure thermal bremsstrahlung 
calculated using the Gaunt factor given in Appendix A. The histogram, which includes the effects 
of Comptonization, shows little or no departure in the Tp = 0.1 cases. The situation, however, 
is different for Tp = 1.0 where the effects of Comptonization are manifest by the exess at high 
energies. The effects of Comptonization and pair annihilation are clearly seen in the high-z 
states. Here the dotted curves represent the sum of annihilation and thermal bremsstrahlung. 
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The Comptonized spectra appear harder over most of the energy range and saturate at energies 
near 30 to form a Wien bump. In addition, the difference between the hard compactnesses of the 
high-z and \ow-z states decreases with increasing and Tp until the two solutions merge as shown 
in Figure 4. This is apparent in Figure 6 by examining the overall levels of the emission of the two 
states. 

In Figure 7 we show the effect of injecting photons with a uniform spatial distribution 
and a blackbody energy distribution at kT = 0.1 keV on the emergent spectra for the high-z 
pair-balanced states of Figure 6. The presence of the soft photon graybody is clearly visible 
and the individual histograms can be identified with values given for Is in each panel according 
the relative levels of the graybody peak. The histograms displaying no soft graybody [Is = 0) 
correspond to the high-z solutions of Figure 6. In each case the soft photon compactness is 
increased to roughly the maximum value permitted by the condition of pair-balance. It is evident 
from this figure that for higher and Tp the plasma becomes increasingly unable to accomodate 
the injection of soft photons. Hence, for ^ 1.0 and Tp ^ 1.0, steady solutions obtained neglecting 
the intrinsic thermal bremsstrahlung of the plasma are not physically realistic. 

In Figure 8 we plot the isotherms in T-lh space, where T represents the 2 and 10 keV spectral 
index here and in the subsequent discussion. Each panel corresponds to a different value of Tp. 
The isotherms are generated by fixing and Tp, varying the value of the parameter Is, and 
determining the values of z for which the plasma is in pair balance. The symbols in this figure 
give the numerically calculated solutions, and the endpoints of the isotherms correspond to the 
Is = solutions shown in Figure 4. The different values of Is at each point can be determined 
from Figure 5. The pair-dominated solutions generally cluster near the high-compactness end of 
the isotherms, and pair-dominated solutions at a given temperature coincide, independent of Tp, 
provided that Tp is not too large. This is because the Comptonization of soft photons and the pair 
balance condition is determined by the Thomson depth r which is ^ Tp when z ^ 1. 

We can divide each panel conceptually into three regimes. In the left-hand portion we find 
the bremsstrahlung regime, where the soft photon energy density is dominated by bremsstrahlung, 
the solutions have few pairs [z <C 1), and the spectra are hard, corresponding to nearly pure 
bremstrahlung solutions. This regime represents the solutions for which Is ^ Iff (see equation 
4-1), which are located near or at the left-most endpoints of each curve. As Is becomes ^ Iff, 
the spectra soften and reach a plateau where the formation of the spectrum is dominated by 
non-bremsstrahlung soft-photon Comptonization and Ih oc Is- This is the linear regime, which 
continues with increasing Is until the Compton opacity due to pairs is comparable to that of the 
ionization electrons. As discussed previously, a second solution with large Compton opacity from 
pairs is found through the condition of pair balance. The solutions in this pair-dominated regime 
are harder than the corresponding pair-dilute solutions for a given value oi Is, and persist even if 
Is = (see Figure 4). 
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Fig. 4. — Hard compactness Ih as a function of for thermal plasmas in pair balance when all 
soft photons are supplied by bremsstrahlung [Ig 0). Curves are labeled by different values of the 
proton optical depth Tp. The dotted and dashed curves respectively show the numerical result of 
Zdziarski (1984) and the analytic result of Svensson (1984) for the case Tp = 1. 
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Fig. 5. — Hard compactness th as a function of for thermal plasmas in pair balance for soft 
photons with compactness tg injected uniformly into a spherical cloud with uniformly distributed 
electrons and pairs. The solutions in the separate panels correspond to different values of the proton 
optical depth Tp. 
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Fig. 6. — Spectra corresponding to the low-z and high-z pair-balanced states of a plasma with the 
parameters© ={0.2, 1.0}, Tp ={0.1, 1.0} and £s=0. The dotted curves show the bremsstrahlung and 
pair annihilation spectra of the plasma without Comptonization. The solid histogram represents 
the fully Comptonized spectra. 
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Fig. 7. — Spectra for the high-z pair-balanced states of a plasma for various values of the soft 
photon injection taken to be a diluted blackbody with kT = 0.1 keV and compactness as indicated 
in the individual panels. The parameters and Tp are the same as those used in Fig. 5. 
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Fig. 8. — Photon 2-10 keV spectral index T as a function of emergent hard compactness Ih for 
thermal plasmas in pair balance. The panels are labeled by the values of proton optical depth 
Tp. The curves represent isotherms with dimensionless temperature given by the symbols in the 
legend. The light shaded region represents the envelope of the entirety of the AGN data in listed 
in Table 1 and the smaller dark region represents the Seyfert 1 "island" discussed in the text. The 
dotted curve delineates the boundary between the pair dominated (z > 1) and pair dilute (z < 1) 
regions, with the pair dominated solutions lying to the left. The vertical dashed line represents the 
efficiency limit of equation (6.4). 
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6. Comparison with Data of AGNs 



We consider the sample of 21 AGNs listed in Table 1, which includes Seyferts, Narrow 
Emission Line Galaxies, and one Broad Line Radio Galaxy. Except for NGC 4151, which we 
classify separately because of its Sy 1.5 to Sy 2 transitions, all 2-10 keV luminosities and spectral 
indices are from the analysis of Ginga data by Nandra & Pounds (1994). They consider the effects 
of photoelectric absorption and reflection by cold matter near the source on the intrinsic source 
spectra. For NGC 4151, we use the luminosity derived by Awaki et al. (1991), who subtracts 
efl"ects of photoelectric absorption by surrounding cold gas on the intrinsic spectrum of NGC 4151. 
Reflection is not important for NGC 4151, as shown by Maisack (1993) and Maisack & Yaqoob 
(1991). We use the measured X-ray variability doubling time scales At quoted by Lightman & 
Zdziarksi (1987) and Done & Fabian (1989) to determine the hard compactnesses Ih through the 
relation 



The middle term in parentheses is a bolometric correction to the 2-10 keV compactness, and 
depends on the spectral index and cutofl" energies Ei and E^- 

Figure 9 shows the location of the 21 AGNs in the V — Ih plane. As can be seen, the intrinsic 
spectral indices cluster in the range 1.5 ^ F ^ 2.5, and the hard compactnesses are found mostly 
in the range 0.01 ^ Ih ^ 10^. The value oi Ih may be a lower bound if the measured At values are 
limited by statistics or incomplete sampling. The error bars indicate the uncertainty in the range 
of Ih due to uncertainty in the values of Ei and E^- The lower limit corresponds to Ei = 1 keV 
and Eu = 50 keV, whereas the upper limit corresponds to to Ei = 10 eV and E^ = 500 keV. Note 
that there is a marked clustering of Seyfert 1 AGNs for 2 ^ ^ 10 and 1.8 ^ F ^ 2.2. We call 
this region the "Seyfert 1 island," but more observations must be made to assess the signiflcance of 
this clustering. (For a discussion of X-ray observations and variability time scales, see Mushotzky 
et al. 1993 and reference therein.) We also see that the NELGs and NGC 4151 evidently have 
smaller values of Ih than the Seyfert Is, although this may be a selection efl"ect caused by the 
larger photoelectric absorption and weaker X-ray luminosity for these sources, which make the 
determination of At more difficult. 

In Figure 8, we plot shaded circles representing the approximate location of the Seyfert 1 
island on the theoretical curves which give F as a function of Ih for thermal plasmas in pair 
balance. We suggest that thermal pair-balanced solutions for Seyfert AGNs are unlikely when 
Tp ^ 0.01 for two reasons. The first is due to the allowed limits on compactness. Clearly not all 
compactnesses are physically possible because, following Lightman & Zdziarski (1987), 



The final inequality holds if the source radiates at the Eddington limit and dissipates this power 
within 10 Schwarzschild radii. Actually the allowed compactnesses can be restricted even further. 




(6-1) 




(6-2) 
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using the efficiency limits of Fabian (1979). Because the radiation escapes from a spherical source 
on the timescale At = R[l + Tp/3)/c, the energy liberated from the matter is bounded above by 
the expression 

LAt ^ —TrR^Up'qmpC^, (6-3) 
3 

where 77 = 10~^77_i is the efficiency of converting the rest mass of matter into radiant energy. 
Using the definition oi I and Tp = npaxR in equation (6-3), one finds that 

4 / 'm„\ T„ T„ , , 

1<^^V(^) = 7707?_i " 6-4 

3 \meJ l + Tp/3 l + Tp/3 

Thus we see that for 10% efficiency, there are no allowed solutions for several high compactness 
AGNs if the thermal plasmas have Tp ^ 0.01. The 10% efficiency limits are shown in Figure 8 by 
the vertical dashed lines. In this argument we ignored the contribution of pairs to the scattering 
opacity in estimating the photon diffusion time At = [1 + Tp/3)R/c because pair densities can 
attain equilibrium on even faster time scales. We point out, however, that the arguments leading 
to equation (6-4) can be weakened by considering other geometries or by extracting energy from 
sources other than the rest mass of the accreting matter (Rees 1982). 

The second reason that the low optical depth solutions are unlikely is due to the observed 
spectral smoothness and relative uniformity of measured spectral indices of Seyfert AGNs, 
compared to results from spectra calculated with small optical depths. When Tp ^ 0.01, the 
allowed solutions in the vicinity of the Seyfert 1 island show distinct spectral structure arising 
from individual scattering profiles, as shown in Figure 3 (see Pozdnyakov et al. 1983; Zdziarski et 
al. 1994). The large variations in spectral index expected in the low optical depth regime have not 
been seen in broadband X-ray measurements of Seyfert AGNs (e.g., Rothschild et al. 1983). 

When Tp ^ 0.1, the values of T become increasingly uniform for a given temperature. We see 
from the location of the Seyfert 1 island that temperatures ~ 0.4 — 0.6 are consistent with the 
data, although the efficiency constraint still requires larger values of Tp in some cases. At still 
larger values of Tp, smaller temperatures are required to give agreement between the observations 
and calculations. Thus we see that if Seyfert AGNs are modelled by thermal plasmas in pair 
balance, then temperatures ^ 0.6, or kT ^ 300 keV, are inconsistent with the data. Moreover, 
these solutions are invariably in the low-pair regime. This is made clear by noting that the AGNs 
lie to the left of the dotted curve indicating z = 1 m Figure 8, which separates the pair-dilute and 
pair-dominated regimes. Thus we predict that if Seyfert AGNs are thermal, then they should show 
spectral cutoffs corresponding to temperatures ^ 300 keV. We also conclude that observations 
of Seyfert AGNs strongly suggest that the plasmas responsible for the high energy emission are 
unlikely to be pair-dominated. The condition of pair dominance was used to reduce the number of 
parameters and establish a one-to-one correspondence between the observables and the physical 
parameters of the source (Ghisellini & Haardt 1994). But our comparison with the inferred 
compactnesses of AGNs shows that AGN data are found in the pair-dilute regime. 

Present observations are not sufficiently sensitive to determine the energy of spectral cutoffs of 
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Table 1. Seyfert Data 



Name 


Type 






log At 


4 


Mrk 335 


Sy 1 


2.39±0.15 


43.5 


3.8 


64±49 


Fairall-9 


Sy 1 


1.94±0.09 


44.5 


5.5 


4.8±2.1 


NGC 3227 


Sy 1 


1.95±0.07 


42.4 


4.3 


2.3±1.0 


Akn 120 


Sy 1 


1.93±0.11 


44.2 


3.4 


270±120 


NGC 3516 


Sy 1 


2.08±0.09 


43.0 


4.0 


4.1±2.1 


NGC 3783 


Sy 1 


2.11±0.07 


43.6 


4.3 


8.2±4.4 


NGC 4051 


Sy 1 


1.84±0.09 


41.5 


3.0 


4.6±2.0 


NGC 4593 


Sy 1 


1.81±0.07 


43.0 


3.4 


3.8±1.6 


MCG-6-30-15 


Sy 1 


2.12±0.04 


43.0 


3.3 


150±82 


NGC 5548 


Sy 1 


1.81±0.02 


43.6 


4.3 


7.8±3.4 


Mkn 841 


Sy 1 


2.30±0.02 


43.7 


5.0 


3.1±2.2 


Mkn 509 


Sy 1 


1.86±0.02 


44.4 


5.7 


2.3±1.0 


NGC 7469 


Sy 1 


1.97±0.05 


43.6 


4.9 


9.0±4.1 


NGC 4151 


Sy 1.5 


1.48±0.04 


42.7 


4.6 


1.2±0.7 


NGC 2110 


NELG 


1.93±0.11 


43.0 


4.9 


0.47±0.21 


NGC 2992 


NELG 


2.30±0.21 


42.7 


5.1 


0.067±0.046 


NGC 526A 


NELG 


1.49±0.21 


43.2 


4.9 


0.027±0.015 


NGC 7314 


NELG 


1.91±0.11 


42.6 


3.3 


9.1±4.0 


NGC 5506 


NELG 


2.08±0.05 


43.2 


5.2 


65±33 


MCG-5-23-16 


NELG 


1.86±0.12 


43.1 


5.0 


0.24±0.10 


3C 382 


BLRG 


1.73±0.11 


44.3 


5.4 


3.4±1.6 



Lx is the 2-10 keV luminosity assuming Hq = 50 km s ^ Mpc ^ and 
Tx is the 2-10 keV underlying power law. 
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Fig. 9. — A plot of the photon spectral index versus hard compactness for the 21 AGN listed in 
Table 1. 
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Seyfert AGNS except in a few cases. OSSE observations of NGC 4151 (Maisack et al. 1993) show 
a strong spectral cutoff at ~ 50 keV which is consistent with our prediction. NGC 4151, however, 
is not a typical Seyfert 1. IC 4329 A is a standard Seyfert 1, but the presence of a reflection 
component makes it possible to flt the OSSE obersations (Madejski et al. 1994) with spectral 
cutofl"s ranging from 300-700 keV. The average spectrum determined from OSSE observations of 
Seyfert 1 AGNs (Johnson et al. 1994) displays a cutofl" at ~ 50 keV, but it is important to take 
into account the reflection component found in such sources (Nandra & Pounds 1994). When this 
is done, analysis of OSSE observations of four radio quiet AGNs, excluding IC4329A, shows a 
spectral softening with average e-folding energy of several hundred keV (Zdziarski et al. 1995). 



We have developed a Monte Carlo thermal Comptonization code which contains accurate 
ep, ee, and e+e" bremsstrahlung and pair annihilation spectra. Pair production through 77 
attenuation is calculated self consistently from the Comptonized angle and energy- dependent 
photon distribution. We also considered internally generated soft photons. We have compared our 
results for the case of pair free plasmas with the analytic expressions for thermal Comptonization 
found in the literature. We constructed a grid of solutions for pair-balanced plasmas and presented 
the results in the spectral index/compactness plane. We then compared our solutions with the 
X-ray data from a sample of 21 Seyfert AGNs and found that pair-dominated solutions are 
inconsistent with the data if the X-ray variability time scale gives an accurate measure of the 
compactness. Furthermore, we predict that the temperatures of central engines of Seyfert galaxies 
are ^ 300 keV. This prediction can be tested with more sensitive gamma-ray observations of 
Seyfert galaxies. 

We thank Dr. A. A. Zdziarski for useful discussions. 



We generate fltting functions accurate to within 8% for the Gaunt factors gi2{x; 0) 
of ep, ee, and e+e" thermal bremsstrahlung using results of numerical integrations of the 
emissivities given by Dermer (1986). The fltting functions reduce to asymptotic forms at 
nonrelativistic temperatures, extreme relativistic temperatures, and in the soft-photon regime. In 
the transrelativistic regime. Stepney & Guilbert (1983) and Haug (1987) also provide tables for ee 
and e+e" thermal bremsstrahlung emissivities, respectively. The form of the Gaunt factors in the 
soft photon regime, a; < 0, approaches the soft-photon Gaunt factors given by Svensson (1984): 



7. Summary 



A. Thermal Bremsstrahlung Gaunt Factor 




(Ala) 
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^+-(-,0) = [-^) e^/eK,{l/Q) 2 ' + ^0 + 20 j e , (Ale) 

where u = a;/0, 77 = 0.56146, ci = 3.42, C2 = 11.2 and C3 = 10.4. Corrections to these expressions 
are made using the functional form 

gr2ix; 0) = gll^{x; 0)/[l + ai2(0) + 6i2(0H. (A2) 

The coefficients ai2(0) and 612(0) are given by 

r 0.130°-^'* < 0.3 
aep(0) = < 0.060 0.3 < < 0.8 (A3a) 

I O.O5O0-°-92 > 0.8 

r 1.6O0°-3i < 0.25 
bep{Q) = < 1.04 0.25 < < 0.5 (A3b) 

I O.880-°-24 > 0.5 

r -O.240°-4O < 0.6 
«ee(0) = < -0.21 0.6 < < 1.5 (A3c) 

I -0.240-°-^^ > 1.5 

, r O.530-°-27 < 1.5 , . 

6ee(0) = < ^ - (A3d) 

I 0.46 > 1.5 

rO.120°-^i < 0.5 
a+_(0) = i 0.082 0.5 < < 1.5 (A3e) 

U. 120-11 > 1.5 

r O.960° i2 < 0.5 
6+_(0) = I 0.680-0-31 0.5 < < 1.5 (A3f) 

I O.6O0-° °95 > 1.5. 

We also give corrections for thermal bremsstrahlung in the hard photon regime, < a; < 100. 
For ep thermal bremsstrahlung with < 0.7 we write the Gaunt factor in the form 

g,p{x; 0) = g:;{x; 0)[1 + a(0) + 6(0)« + c{e)u% (A4) 

where the nonrelativistic Gaunt factor is 

g:^{x;Q) = —exp{x/2Q)Ko{x/2Q), (A5) 

TT 

(Green 1959) and the coefficients a(0), 6(0) and c(0) are given by 

a(0) = 1.60i-^ (A6a) 



For > 0.7 we have 



where 



KQ) _ fO.360°-92 < 0.4 

1 0.14 O.4<0<O.7 

c(0) = 0.530^-^^. 

gepix; 0) = 1 + a(0) + 6(0)w + c{e)u\ 



2.30 



0.71 



{ 



0.2203-1 O.7<0<1.O 

O.2201-® 1.0 < < 2.0 

0.3500.98 > 2.0 

0.0990^ '* 0.7 < < 2.5 

0.14011 > 2.5. 



a(0) = 

6(0) = 
c(0) = 

For ee thermal bremsstrahlung 

gee{x; 0) 

where 

a(0) = 
6(0) = 

c(0) = 

In the case of e+e" thermal bremsstrahlung with < 0.45 we have 



+ a(0) + 


6(0)u+ c(0)u2 


3.001-1^ 


< 1.5 


33300.814 


> 1.5 


0.4101-6 


< 0.5 


0.2801-1^ 


> 0.5 


0.4602-4 


< 0.6 


0.3401-"^^ 


0.6 < < 1.5 


0.450106 


> 5.5 



g+-{x; 0) 



0.9194 

with gl''_{x; 0) = 2V2gl^^{x; 0) (Haug 1985) and 



gl'_{x; 0)[1 + a(0) + b{Q)u + c{Q)v 



a(0) 



' 1O.O102-'19 < 0.15 
4.702-°9 0.15 < < 0.3 
2.0301-36 0.3 < < 0.45 
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6(0) = -O.2240°-^3 (A12b) 

c(0) = 1.290^-® (A12c) 
And finally, for > 0.45 we have 

g+-{x; 0) = 1 + a(0) + b{Q)u + c{Q)u'^, (A13) 

where 

r 7.30°-SS 0.45 < < 1.5 /Ai. X 

6(0) = 0.410^-3^ (A14b) 

c(0) = 1 0-^40-^ O.45<0<2.O (^^4^) 

^ ^ ll.O80°-98 0>2.O ^ ^ 
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